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The Born-Oppenheimer electronic wavefunction ^R^(r) picks up a topological phase factor ±1, 
a special case of Berry phase, when it is transported around a conical intersection of two adiabatic 
potential energy surfaces in i^-space. We show that this topological quantity reverts to a geometric 
quantity if the geometric phase 7 = ^ Im($i?| is evaluated with the conditional 

electronic wavefunction ^R{r) from the exact electron-nuclear factorization ^R{r)x{R) instead of 
the adiabatic function A model of a pseudorotating molecule, also applicable to dynamical 

Jahn-Teller ions in bulk crystals, provides the first examples of induced vector potentials and molec¬ 
ular geometric phase from the exact factorization. The induced vector potential gives a contribution 
to the circulating nuclear current which cannot be removed by a gauge transformation. The exact 
potential energy surface is calculated and found to contain a term depending on the Fubini-Study 
metric for the conditional electronic wavefunction. 
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I. INTRODUCTION 


The Born-Oppenheimer approximation underlies most 
calculations in condensed matter physics and chemistry. 
Examples include thermal conductivity, lattice-mediated 
relaxation of excited electrons and optical properties of 
materials, as well as molecular scattering and rovibronic 
spectroscopy. Since nuclei are much heavier than elec¬ 
trons, one can get a good approximation to the electron- 
nuclear wavefunction ^(r, R) ^ by as¬ 

suming the nuclei are frozen and solving an electronic 
Schrodinger equation with the R-dependent Hamiltonian 
fjBO — ifin, which is the full electron- 

nuclear Hamiltonian with the nuclear kinetic energy re¬ 
moved. The eigenvalue defines an adiabatic potential en¬ 
ergy surface 

+ Fee + Gn + (1) 

which is then used in the nuclear Schrodinger equation 
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+ x^^{R) = Ex^^iR). ( 2 ) 


the full wavefunction must be a single¬ 
valued function of R, multivaluedness of implies 

that ^R^{r) is also multivalued, so it does not return 
to its original value if transported along a closed path 
in R-space encircling a conical intersection, but instead 
changes sign. This sign change is due to the Longuet- 
Higgins phase.^ It is a special case of the Berry phase^ 
because it only takes the values 0 and tt. 

Multiplication by a Dirac phase factor exp ^ f • dR^ 
which compensates the sign change makes ^R(r) single¬ 
valued. If this choice of phase is made, lm{^ r\\/r) is 
no longer zero, and Eq. (2) must be replaced by^’^ 
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where is the induced vector 

potential introduced by Mead and Truhlar.^ 

Nonadiabatic coupling between electronic eigenstates 
causes corrections to the Born-Oppenheimer approxima¬ 
tion, which are usually included through the expansion^”^ 


<i!{r,R)=Y,^n{r,R)Xn{R), (4) 
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The adiabatic potential energy surface is an extremely 
useful concept which implicitly encapsulates all electronic 
terms (kinetic Tg, interaction I4e and electron-nuclear 
coupling Ven) and the nuclear interaction Vnn ia a single 
scalar function E^^{R) under the assumption that the 
electronic wavefunction <I>^^(r) stays in the ground state 
of the electronic Hamiltonian for all values of the nuclear 
coordinates R; we use the notations R = (Ri,R 2 ,...) 
and r = (ri,r 2 ,...). 

A curious feature of the Born-Oppenheimer approxi¬ 
mation is the occurrence of conical intersections between 
the potential energy surfaces of two or more electronic 
eigenstates in some polyatomic molecules.The factor 
{R) is then multivalued due to the nonanalyticity of 
the potential £^^{R) at the point of intersection. Since 


i.e. through a sum of Born-Oppenheimer-like factors, one 
for each stationary state <I>n(r, R). Surprisingly, it is not 
actually necessary to depart from the Born-Oppenheimer 
single-product form to include nonadiabatic effects. In 
fact, the exact wavefunction can be factored into a sin¬ 
gle product ^R{r)x{R) called the exact electron-nuclear 
factorization}^ Like the Born-Oppenheimer ansatz, 
the exact factorization can be used to define scalar and 
vector potentials £{R) and A^{R). 

Using these exact potentials in place of 
in Eq. (3) yields the exact nuclear density and nuclear 
current density of the state ^(r, Since it integrates 
all nonadiabatic electronic effects into a single potential 
energy surface and vector potential, the exact factoriza¬ 
tion provides an intuitive and economical description of 
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quantum nuclear dynamics. From this standpoint, the 
following two questions are relevant. How does the exact 
potential energy surface differ from the adiabatic ones, 
e.g. do conical intersections persist? Do the vector po¬ 
tential and Longuet-Higgins phase connected with the 
nonanalyticity at the point of conical intersection remain 
nonzero in the exact factorization? If the vector poten¬ 
tial can be made to vanish by a gauge transformation, 
then the potential energy surface is all one needs to de¬ 
scribe the nuclear motion. If instead the molecular geo¬ 
metric phase is nontrivial, then the nuclear Schrodinger 
equation must contain induced vector potentials. Recent 
work found a case in which the Longuet-Higgins phase of 
TT becomes zero in the exact factorization.^^ However, the 
model studied does not have the degeneracy of the classic 
Jahn-Teller models^^ of pseudorotating molecules^’^’^^”^^ 
and transition metal ions in bulk crystals^^”^^ and there¬ 
fore leaves room for further investigation. 

Within the Born-Oppenheimer approximation. Berry 
phase effects have been extensively studied in pseudoro¬ 
tating molecules such as hydrogen exchange 

reactions^^~^^ and fullerene ions and crystals.^^”^^ Berry 
phase effects are also relevant to some dynamical Jahn- 
Teller systems which have been investigated recently. 

In many of these cases, the exact factorization could 
provide an interesting alternative to conventional ap¬ 
proaches based on the adiabatic expansion of Eq. (4). 
Because the nuclear wavefunction is determined by a sin¬ 
gle potential energy surface and vector potential, rather 
than an infinite set of coupled adiabatic potential en¬ 
ergy surfaces, the exact factorization scheme might prove 
to more efficient than traditional approaches to cou¬ 
pled electron-nuclear dynamics if accurate approxima¬ 
tions can be found for £ and A^. 

A two-mode vibronic model has been studied in the 
context of the exact factorization to explore features of 
the exact potential energy surface, such as spikes that oc¬ 
cur near nodes of the adiabatic nuclear eigenfunctions.^^ 
Such spikes were previously observed in the ID poten¬ 
tial energy surfaces of diatomic molecules.Here, we 
focus on the exact vector potential in a model pseudoro¬ 
tating triatomic molecule known to have a degenerate 
ground state and nonzero Longuet-Higgins phase in the 
adiabatic approximation.^^ The presence of degeneracy, 
which is maintained in the exact approach, is pivotal: 
the exact factorization can be applied to any state in the 
ground state manifold, and the potential energy surface 
and vector potential will depend on which state is chosen. 
Any choice breaks the unitary symmetry of the ground 
state manifold and leads to a corresponding symmetry 
breaking of the exact potential energy surface and/or 
vector potential. Our main result is that the discrete 
topological Longuet-Higgins phase, 0 or tt, becomes a 
full geometric phase in the exact approach. 

Section H introduces our model pseudorotating tri¬ 
atomic molecule. After solving the full electron-nuclear 
Schrodinger equation in Sec. HI, we evaluate the vector 
potential and molecular geometric phase in Sec. IV, and 



FIG. 1: Two sets of unit vectors an and bn used to specify 
the geometry of the distorted triatomic molecule. 

the exact potential energy surface in Sec. V. The Berry 
curvature and a Riemannian metric derived from the con¬ 
ditional electronic wavefunction are unified in a quantum 
geometric tensor in Sec. VI. 

II. MODEL PSEUDOROTATING MOLECULE 

Our model pseudorotating molecule consists of three 
hydrogen-like atoms, whose nuclei are assumed to be 
distinguishable to avoid the complication of nuclear ex¬ 
change symmetry. This models a molecule like LiNaK, 
but we further assume the nuclear masses are equal. We 
will follow, as closely as possible, the notations of Ref. 34, 
where the model was introduced. The well known E ^ e 
Jahn-Teller system can be obtained by the truncation to 
two electronic levels, which would be a good approxima¬ 
tion when the electronic level spacing is much larger than 
the characteristic vibrational energy. 

The electronic degrees of freedom are described within 
the truncated Hilbert space formed by three valence elec¬ 
trons occupying three 5-like orbitals, one per atom, while 
the full real space i?-dependence of the nuclear states is 
retained. The Hamiltonian is 

^ = + (5) 

where is the nuclear kinetic energy, Wn is the nuclear 
interaction energy and 

Hen — ^ ^ H.C.^ (6) 

ncr 

subsumes electronic kinetic energy and electron-nuclear 
coupling. The electron-nuclear coupling is represented 
in the i?-dependence of the hopping amplitudes tn,n+i- 
Electron-electron interactions are neglected; their effect 
has been studied in the adiabatic limit. 

The key simplifying assumption of the model is the 
truncation of the Taylor expansions of Hen and Vnn with 
respect to distortions of the equilateral triangle geometry 
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(see Fig. 1). For the hopping amplitude tn,n+i(^) in 
Eq. ( 6 ), we keep only the linear term 


— to — -^(|Rn+l “ Rn | ~ VSRq) 

— to —(U72+I U72) • (S-n+l 


where g is the electron-nuclear coupling constant, VSRo 
is the internuclear separation in the equilateral config¬ 
uration, Un is the deviation of nuclear coordinate R^ 
from its equilateral position, i.e. = R^ — and 

a.n is the set of unit vectors shown in Fig. 1 . We choose 
to = 1 to be our characteristic unit of energy. The term 
Vnn models internuclear repulsion and some part of the 
electron-nuclear attraction. We assume harmonic spring 
interactions by expanding to second order in U^: 


Vnn — (|Rn+l — Rn | ~ \f^Ro) 

n=l 

= ^Ql + ’^Q^ (8) 

where Qi and Q are internal nuclear coordinates defined 
below, and we have introduced the effective spring con¬ 
stants Ki = 3K and K = |i^. 

The three nuclear coordinates are determined by three 
center of mass coordinates, three Euler angles and three 
internal coordinates. The nuclear coordinates in the lab¬ 
oratory frame are 

R lab T3 I 'T^'Dmol /q\ 

where Rem is the nuclear center of mass, T is an SO(3) 
rotation and R™°^ are molecular frame coordinates. The 
molecular frame coordinates can be parameterized by 
three normal mode amplitudes (Qi,Q 2 ,Q 3 ); Qi is the 
breathing mode; Q 2 and Q 3 are the conventional sym¬ 
metric and asym metric be nding modes.In polar co¬ 
ordinates Q = r] = tan“^(( 53 /Q 2 ), we have 


Rr' = + Qi)a„ + QK 


= cos 


27rn\ . /27rn\ 


b„ = ( cos (r? - j, sin (r? 


The 7 ^-dependent unit vectors are shown in Eig. 1. 
The nuclear kinetic energy operator is 



2M\QdQ\^dQ)^ Q^dT]^^ dQl)' ^ 


The total nuclear mass M. = 3M appears here because 
(Qi, Q, T]) describe collective nuclear motion. As one fur¬ 
ther simplification, we neglect the last term in Eq. (11) 
and assume that Qi is frozen to its equilibrium value. 
Since Qi is a fully symmetric mode, this simplification 
will have no effect on our qualitative conclusions concern¬ 
ing the molecular geometric phase. The Hamiltonian is 
now fully defined, and we proceed to calculate its eigen¬ 
states. 


III. EXACT SOLUTION OF THE MODEL 

The model will be solved by exact diagonalization af¬ 
ter the Hamiltonian matrix elements are calculated in an 
electron-nuclear product basis. 

We start by defining a complete set of electronic states. 
The truncated 3-electron Hilbert space is 20-dimensional, 
but since and Sz commute with the Hamiltonian, 
we focus on the 8-dimensional sector with spin quan¬ 
tum numbers S = ^ and 5'^ = The basis states 
are constructed from the ket |core) representing the in¬ 
ert core electrons by acting with the creation operators 
of single-particle orbitals |0/c) = ^ X]n=o 1 2 

k = 1,0,—1, e.g. |t1t 0) = where the 

position of the spin in the ket represents the orbital it oc¬ 
cupies with the ordering convention k = 1,0, —1. Thus, 
we choose the following eight basis states: 

|a) = |tt0) 

l^) = |ott) 

k) = |tot) 

|^) = |t Of) 

k) = |ott) 

I/) = It to) 

i<?) = -te t; t) -111;) - u 11)) 

l/i) = -T(^oitit) + ittb-Utt)). (12) 

The nuclear basis states in the polar coordinate repre¬ 
sentation are with radial and azimuthal 

quantum numbers n and m. As a complete basis of radial 
functions, we choose the normalized radial eigenfunctions 
of the isotropic 2D harmonic oscillator 

^ |m|)! 

where N = {n — |m|)/2, A = and L^m is the 

associated Laguerre polynomial, defined by the differen¬ 
tial equation 

xy" -h (m + 1 — x)y' + Ny = 0. 

The next step is to evaluate the Hamiltonian matrix 
elements in the electron-nuclear basis |anm) = |a)|nm): 

(ainimi|^|a2n2m2) = ^aias + Kn|^2^2) 

+ (ainimi|^en|«2^2^2). (13) 

All matrix elements can be evaluated analytically and are 
reported in App. A. 

Solving the model by exact diagonalization confirms 
that the ground state is doubly degenerate for certain 
values of the parameters M, K and g. Before going on 
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to investigate the molecular geometric phase, it is in¬ 
structive to examine the symmetries of the model. 

Threefold symmetry is responsible for several special 
properties of the model. A combined symmetry operator 
C 3 = where Cse is a threefold permutation of 

the electrons (0 ^ 1 ^ 2 ^ 0) and is a threefold 
pseudoTotaiion of the nuclei (77 ^ 77 — 27r/3), commutes 
with the model Hamiltonian in Eq. (5). This symmetry 
is what remains of the rotational symmetry of the origi¬ 
nal real-space Hamiltonian in our model in which overall 
molecular rotations are not included. To see this, we note 
that if Cs is combined with a three-fold permutation of 
the nuclei 0 ^ 1 —> 2 ^ 0, denoted by Csn, then the 
combined operation C^C^n corresponds to a 27r/3 rota¬ 
tion of the entire molecule, which is subgroup of the full 
rotational symmetry of the original Hamiltonian. 

The operator Cse shifts the electrons forward by one 
site, e.g. CseltItO) = | 0t1t) foi* arbitrary ket expressed 
in the (0,1,2) site basis. The basis functions in Eq. (12) 
were chosen to be eigenstates of C^e'- |c), |e) have 

eigenvalue |6), jd), |/) have eigenvalue 

1^) and \h) have eigenvalue 1. Similarly, the operator Csr] 
rotates the nuclei by —27tI3 in (Q 25 Q 3 ) space as shown 
in Eq. (16) below. 

The threefold symmetry group generated by Cs has 
two irreducible representations - a symmetric singlet A 
and a doublet E. Being a doublet, the ground state be¬ 
longs to the E representation. We can choose two orthog¬ 
onal states from the ground state manifold that trans¬ 
form into themselves up to a phase under Cs. These 
states have eigenvalues and will be labeled |T±) 

with +/— indicating counterclockwise/clockwise nuclear 
current. Only \anm) states with the same eigenvalue of 
Cs are coupled by the Hamiltonian. Hence, the state 
|T+), for instance, has the following structure: 

|T+) = r,|a) + E^r,\b) + Tele) + e^^T^ld) 

+ Tele) + + Tf,\h)), (14) 

where = Ta{Q,T]) are periodic functions of 77 with 
period 27r/3, i.e. 

CX) 

Ta{Q,V)= E ra,3m(Q)e'='™’'. (15) 

m= — oo 

The state |T_) is the complex conjugate of |T+). The 
operator acts on the nuclear functions as 

C3n^aiQ,v) = r«(Q,C3y7?) = r«(Q,r/ + 27r/3). (16) 

Thus, we immediately verify that CsjT-i-) = 

To see how double-valued electronic and nuclear wave- 
functions, and x(i?), can emerge from the single¬ 

valued function in Eq. (14), consider the simultaneous 
M ^ oo and ^ ^ 0 limit with the condition hCl/A = 
const^ where A = ^^/2/C is the Jahn-Teller stabilization 


energy. In this limit, it can be shown that 

= p(g)e*>'/2(^e-*’'/2|a) + (17) 

The function in parentheses is the real-valued electronic 
Born-Oppenheimer function |4>^^), showing the charac¬ 
teristic sign change when 77 increases by 27r. The nuclear 
factor is xiQ^v) = where p{Q) is approxi¬ 

mately equal to a harmonic oscillator wavefunction cen¬ 
tered at the minimum of the potential — gQ. The 
lowest order corrections to Eq. (17) in the M ^ 00 limit 
will be investigated in further detail elsewhere. 

Instead of |T+) and |T_), one can choose two real¬ 
valued orthogonal states from the ground state manifold, 
e.g. |T^) = \^Re|T+) and |T^) = A/ 2 Im|T+), character¬ 
ized by their parity g/u under the reflection Qs ^ — Qs- 
In the above M ^ 00 , ^ ^ 0 limit, these reduce to 

1^,) ^ ^(1 + COST?)lot 0igt) + ^ sinr? lot Ot u^) 

l^'u) -S' ^^^(1 -cos?7)|0t04,ut) + d^sin7?|0t04,5^), 

where | 0 t 0 ^ ^t) = with the following single 

particle orbitals: 

|0o) = |0/c=o) 

\(pg) = V2Re\4>k=+i) 

|(/.„) = V21m\<Pk=+i)- (18) 

|T^) and are nonadiabatic electron-nuclear counter¬ 
parts of the electronic states | 0 t 0 | g^) and | 0 t 0 | u^). 
Our numerical results are consistent with those obtained 
in the two-level approximation with the electronic state 
constrained to the {| 0 t ^), | 0 t subspace. 

Symmetry analysis tells us that all of the eigenstates 
transform as either A or E irreducible representations, 
but only the Hamiltonian matrix elements can determine 
their energetic ordering and whether the ground state has 
A or E symmetry. In the following, we choose parame¬ 
ters such that the ground state is degenerate {E symme¬ 
try). Two features of the model are chiefly responsible for 
the ground state degeneracy. Eirst, fermionic symmetry 
constrains the third electron to occupy the degenerate 
k = El orbitals after the lowest energy k = D orbital is 
occupied twice. Second, electron-nuclear coupling lowers 
the energy of the pseudorotating E states with respect 
to A states. 

The nuclear probability density is 
|x(i?)|2 = (T|T), 

= ^|r„(i?)|2, (19) 

a 

where (• • • )r denotes the inner product on the electronic 
Hilbert space only. In Eig. 2, \x{R)\‘^ is shown as a 
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pseudorot at ing molecule. It is instructive to review the 
derivation^^ of the geometric phase for an arbitrary open 
curve R{s), parameterized by s G [si, 52 ], from Ri to R 2 . 
One starts from the Pancharatnam phase Arg(^i^^ \ ^R 2 )^ 
which defines a unique relative phase between the end¬ 
points, provided \^r^) and |^i? 2 ) orthogonal.^^ 

The open-path geometric phase is then conventionally de¬ 
fined to be the remainder after subtracting the dynam- 
ical phase, {^R:^s)\H{s)\^R(^s))ds, from the Pan- 

charatnam phase.Here, H{t) is an auxiliary Hamilto¬ 
nian that drives the electronic wavefunction along the 
path |^i?(s)) in Hilbert space, i.e. |^i?(s)) is the solution 
of the Schrodinger equation ihds\^) = H{s)\^). Since 
our definition = h lm{^ r\\/r) differs by a sign 
from the conventional definition = ih{^ r\\/r), we 
have introduced a minus sign in the following definition 
of the exact molecular geometric phase: 


FIG. 2: Nuclear probability density |x|^ for (a) the current- 
carrying state |^+) and (b) the even parity state |^^). Lower 
panels (c,d) show the corresponding potential energy surfaces. 
Peaks in \x\^ in (a) correspond to obtuse triangle geometries. 
Model parameters are M = 24, A = 0.6 and g — 1.20 and 
energy is measured in units of the hopping parameter to- 

function of (Q 2 ,Q 3 ) for model parameters M = 24, 
A = 0.6 and g = 1.2. The nuclear wavefunction was ex¬ 
panded over the basis {Qri\nm) = with 

n = 0,1,... 13, m = —n, —n + 2,..., n and \m\ < n. 
The peaks at 77 = (tt/S, tt, 57 r/ 3 ) correspond to the three 
distinct obtuse triangle configurations, which confirms 
that our model is qualitatively consistent with adiabatic 
results for The E 0 e Jahn-Teller model 

with linear electron-nuclear coupling oc gQ^ c.f. Eq. (A 2 ), 
has cylindrical symmetry with respect to rotations in the 
(Q 25 Q 3 ) plane. Quadratic and higher-order terms break 
cylindrical symmetry. In our model, the coupling of \a) 
and I b) to the other six electronic states lowers the cylin¬ 
drical symmetry to Cs even though we have only linear 
electron-nuclear coupling. The potential energy surfaces 
and electron-nuclear coupling of real molecules generally 
have significant anharmonic contributions,^^ which we 
are neglecting. 

Given |xp, the marginal nuclear probability amplitude 
can be written as 

X{R) = (20) 


1 

7 =- - J 

= -Arg($iji + ^ / (21) 

^ jRi 

with an implicit sum over nuclei ja. Equation (21) extends 
the familiar geometric phase to open paths in a way that 
maintains gauge invariance under i?-dependent gauge 
transformations of \^r). Like the Aharonov-Anandan 
phase,it does not assume an adiabatic approximation. 

Equation (21) is also valid in the Born-Oppenheimer 
approximation. In that approximation, there are two 
choices of gauge for which the Longuet-Higgins phase tt 
can be simply understood.Eirst, if the phase of ^i?(r) 
is chosen such that ^R{r) is real-valued for all i?, then 
A^ = 0 and the second term of Eq. (21) vanishes. Hence, 
the phase tt comes from the first term and the sign change 
of ^i?(r), i.e. ^R^{r) = —^R^{r) (multivaluedness). Al¬ 
ternatively, if ^i?(r) is chosen to be single-valued, then 
the first term vanishes for Ri = R 2 and the second term 
gives the phase tt because there is no gauge for which 
the vector potential is zero everywhere along the path. 
Applying i?-dependent gauge transformations to ^i?(r) 
and x{R) changes the first and second terms of Eq. (21), 
but their sum remains the same. 

Using the definition of the conditional electronic wave- 
function ^i?(r) = T(r, R)/x{R)i the vector potential can 
be expressed as 


where S{R) is an arbitrary function of R. Changing the 
phase according to S{R) S{R) + A(i?) implies a gauge 
transformation A^ ^ A^ — V^A of the vector potential. 


IV. MOLECULAR GEOMETRIC PHASE 

With the exact solution in hand, we can demonstrate 
that the geometric phase of ^r is nonzero for our model 


, IxP 

The first term of Eq. (22) is gauge invariant. It is respon¬ 
sible for the geometric phase when the exact factorization 
is applied to our model system. 

Since all three nuclear coordinates R^ are determined 
by (Q, 77 ), it is convenient to define the vectors R = Qeg 
and A = AqQq + in order to simplify the notation. 
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FIG. 3: Nuclear current vector field (arrows) superimposed 
on a color scale plot of the Berry curvature Bq^q^ for |^+). 
The green circle shows the path on which the geometric phase 
(Fig. 4a) and various cross-sections (Fig. 6) are calculated. 


For a closed path C, the geometric phase is 

7 =I^A-dR=I^ [AgdQ + A^Qclri) , (23) 

where Aq = 

Ground state degeneracy has important consequences 
for the molecular geometric phase and vector potential. 
Using the parameters ( 6 >, ip) to express an arbitrary state 
in the degenerate ground state manifold, we have 

1^) = cos + sin (24) 

When ^ = 0, the state has positive (counterclockwise) 
nuclear current in the (Q 2 , Qs) plane as shown in Fig. 3. 
For 6 > = TT, the current is equal and opposite. Thus, 0 
tunes the current continuously between its maximum val¬ 
ues at the poles of the Bloch sphere.Figure 4a shows 
the geometric phase calculated along the circular path 
{Q = 0.54, T] = [0,27r]} as a function of 0. Unlike the 
Longuet-Higgins phase, the exact geometric phase is not 
quantized to 0 or tt and varies in proportion to the nu¬ 
clear current carried by the state. The exact factorization 
scheme can also be applied to adiabatic electron-nuclear 
eigenstates such as the limiting function in Eq. (17). In 
this case as well, the geometric phase can be tuned by 
forming superpositions of degenerate states analogously 
to Eq. (24). Nevertheless, the geometric phase of such a 
state, 7 ^^ = TTCOs^, is still a topological quantity, which 
can only take the values with integer n. 

For 0 = 7 r/ 2 , |^) is real and the current and geometric 
phase vanish for any value of p, since it is always possible 


to choose a gauge such that A is identically zero. How¬ 
ever, C 3 symmetry is broken and the nuclear probability 
density is displaced in the f] = p direction of the (Q 25 Qs) 
plane. For example, p = 0 gives = v^Re|^+) for 
which |x(i?)p is offset in the 7 = 0 direction; cf. Fig. 2b. 
The value p = ir gives the odd state |^u)- 

The topological Longuet-Higgins phase depends only 
on whether the path C encloses a conical intersection, 
otherwise it is path independent. This is not the case for 
the exact molecular geometric phase. Figure 4b shows 
the value of the geometric phase calculated on a circular 
path as a function of the radius Q. The geometric phase 
is path dependent because via Stokes’ theorem it depends 
on the net flux of the effective magnetic field B = V x A 
through the surface S bounded by C: 

= Adq‘' (25) 

where the antisymmetric wedge product is used to write a 
general expression in terms of arbitrary coordinates 
The Berry curvature 

= hlm{d^^R\dj,^R) (26) 

is related to the field strength by Bq^ = with 

= Q 2 and q^ = Q 3 . Since B^j, is an antisymmetric ten¬ 
sor it has only one independent element Bq^q^ which is 
plotted in Fig. 3. For |^+) in the M ^ 00 limit, the mag¬ 
netic flux becomes localized at the origin in the (Q 2 , Qs) 
plane and approaches ^S{Q 2 )S{Q 3 )^ thereby recovering 
the well known adiabatic result. In the three-dimensional 
nuclear coordinate space (Qi, Q 2 ^Q 3 ), the flux coincides 
with the line (generally, a submanifold of codimension 2 ), 
parametrized by Qi, on which the two adiabatic poten¬ 
tial energy surfaces undergo a conical intersection. The 
magnetic field is thus equivalent to that of an infinitesi¬ 
mal flux tube carrying flux h/2 and its Aharonov-Bohm 
phase shift is responsible for the discrete Longuet-Higgins 
phase TT. In the (Q 25 Q 3 ) plane, the vector potential of 
the flux tube is equivalent to that of a magnetic monopole 
with charge ^ = 1/2 at the origin. 

In the exact factorization, the flux tube (or monopole) 
gets smeared out over a finite area, i.e. the point-like flux 



FIG. 4: Molecular geometric phase (a) calculated along the 
circular path in Fig. 3 versus the mixing angle 0 in Eq. (24) 
(b) as a function of path radius Q for ^ = 0. 
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QhWHM (^) QhWHM (b) 




FIG. 5: Dependence of the characteristic radius Qhwhm on 
M and g. Calculations were performed for (a) K — and 
g — 0.01 and (b) K — 2/9 and M = 10^°. 


is replaced by an extended flux density l3{Qi,Q2,Q3) 
that satisfies JJ B{Qi^Q2^Q3)dQ2dQs = h/ 2 . Unlike 
the adiabatic case, the geometric phase vanishes if the 
path encircling the origin is shrunk to a point. On the 
other hand, if the radius of the path is taken to infinity, 
the geometric phase approaches tt, as shown in Fig. 4b, 
since then all of the flux is enclosed. Similarly, the geo¬ 
metric phase calculated on a path with any finite radius 
Q approaches tt in the limit M ^ oo, since the character¬ 
istic radius of the flux tube tends to 0. The characteristic 
radius can be defined as the half width at half maximum 
of the peak in Bq^q^ at the origin. Figure 5a shows that 
Qhwhm decreases as as M ^ oo, consistent with 

the result Qhwhm ^ ^^1/2 of an asymptotic analysis 
to be reported elsewhere. Figure 5b shows that Qhwhm 
varies as if M is sufficiently large. 

For a complex Hamiltonian with three-dimensional 
slow parameter space (X, F, Z), such as Berry’s original 
two-level example,^ degeneracies of the adiabatic poten¬ 
tial energy surfaces occur at isolated points and the geo¬ 
metric phase can be calculated as the flux of a magnetic 
monopole with charge ^ = 1/2 located at the point of 
degeneracy. If the total Hamiltonian includes a kinetic 
energy operator for the slow variables and the separa¬ 
tion of fast and slow variables is made using the exact 
factorization instead of the Born-Oppenheimer approx¬ 
imation, then the geometric phase can be calculated in 
the same way but the point monopole will be replaced by 
a smeared-out magnetic charge density p(X, F, Z) which 
integrates to 1 / 2 . 

The nuclear current density of the state |T) is 

J = M-i/^Im(T|VT)^ 

= 7W-i(ftImx*Vx + A|xH. (27) 

For the state |T+), the nuclear current circulates in the 
positive direction around the origin of the (Q 2 , Q 3 ) plane, 
see Fig. 3, while the electronic particle current (not to be 
confused with the charge current) circulates around the 
three-site ring in the positive direction 0 ^ 1 ^ 2 ^ 0 . 
In our model, the electronic circulating current is 


Looking at the spin-resolved electronic current 

Je<j — ^ ^ H.C.^ (^ 9 ) 

n 

reveals that spin up electrons circulate in the positive 
direction, while somewhat fewer spin down electrons cir¬ 
culate in the negative direction. 

To summarize, when the Born-Oppenheimer product 
is replaced by the exact factorization, the molecular ge¬ 
ometric phase is no longer quantized and can take any 
value between 0 and 27r even though the Hamiltonian is 
real valued. Nonzero geometric phase implies that 
cannot be gauged away and hence gives a nontrivial con¬ 
tribution to the nuclear current density. The curl of 
is an induced magnetic field, which can be viewed as 
the field a smeared-out Aharonov-Bohm flux tube. Spin- 
orbit interactions cause a similar spreading in the adia¬ 
batic case.^^ However, the smearing effect we have stud¬ 
ied here is different because it is a nonadiabatic effect 
of the nuclear kinetic energy operator that occurs even 
when spin-orbit interactions are neglected. 

Since our model takes into account the full eight¬ 
dimensional basis of three-electron states as opposed to 
only the two lowest energy states of the E^e Jahn-Teller 
model, electron-electron interactions can be represented 
in the model Hamiltonian. We have performed calcula¬ 
tions with a Hubbard term U 

and verified that it does not change our main results. 

We have neglected nuclear exchange symmetry and ro¬ 
tations, and therefore the question arises whether the 
molecular geometric phase will survive if they are in¬ 
cluded in the problem. Since the occurrence of nontrivial 
molecular geometric phase depends critically on the de¬ 
generacy of the state, we end this section by commenting 
on the implications of some known degeneracies of the 
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FIG. 6: One-dimensional circular slice (Q = 0.54, g = [0, 27r]) 
of (a) Ixl^ in Fig. 2a and (b) S{R) in Fig. 2c. 
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FIG. 7: Born-Oppenheimer energy surface 
(red) and the potential term {^r\H^^\^r) (blue) for the 
state |^+). The exact potential energy surface is the sum of 
{^r\H^^\^r) and the metric term in Fig. 8. 

general electron-nuclear problem for the molecular geo¬ 
metric phase and induced vector potential. The energy 
eigenstates of any system of electrons and nuclei can be 
chosen to be simultaneous eigenstates of the total angular 
momentum operators Jz^ due to the isotropy of space. 
The resulting 2 J + 1 degeneracies are vital for the exact 
molecular geometric phase because they make it possible 
to construct current-carrying eigenstates from complex 
superpositions of degenerate states. If the molecular ge¬ 
ometric phase of a current-carrying state is nonzero, the 
nuclear Schrodinger equation must contain induced vec¬ 
tor potentials which contribute to the nuclear current. 
Whether the molecular geometric phase is appreciably 
different from its adiabatic value in a particular system 
depends on a number of factors including the strength 
of the electron-nuclear coupling (represented by g in our 
model) and the rigidity of the molecule. Systems in which 
the electronic state is nonadiabatically excited by large 
amplitude nuclear motions and floppy molecules, such 
as pseudorot at ing molecules or highly excited molecules, 
are more likely to have a molecular geometric phase that 
differs significantly from the adiabatic value. 


V. POTENTIAL ENERGY SURFACE 

The Schrodinger equation for the nuclei in the exact 
factorization has the same form as Eq. (3) except the adi¬ 
abatic potential energy surface and vector potential 
A^^ are replaced by their exact counterparts. The exact 
potential energy surface can be expressed as^^ 

£{R) = 

R ^ R ^ 


where = H — Equation 30 is analogous to 

the expression for the adiabatic potential energy surface 
including the Born-Huang nuclear gradient terms. 

The potential energy surface for the nuclear factor 
of the current-carrying state |^+) is plotted in Eig. 2c. It 
appears to have rotational symmetry with respect to the 
pseudorotational angle 77 , however the one-dimensional 
cut along the circle Q = 0.54 in Eig. 6 reveals a weak 
threefold symmetric warping consistent with |x+p. 

The potential energy surface for the nuclear factor Xg 
in Eig. 2d has a high and narrow barrier, which is cut off 
by our choice of scale. As mentioned in the introduction, 
this barrier correlates with the Cs symmetry breaking of 
14^^). The potential energy surface shows such a barrier 
for any real-valued |4^), and its direction is determined by 
the angle (p in Eq. (24). Similar to what was found for di¬ 
atomic molecules^^”^^ and a two-mode vibronic model, 
the barrier in £{R) corresponds to nuclear configurations 
where |x(i7)p drops close to zero. 

The exact potential energy surface is compared with 
the Born-Oppenheimer surface in Eig. 7. The cusp of the 
Born-Oppenheimer surface due to the conical intersection 
at the origin gets smoothed out in the exact surface. 

VI. QUANTUM GEOMETRIC TENSOR 

The last two terms of Eq. (30) can be combined into a 
single quantity 

^geo = (31) 

where is the inverse inertia tensor appearing in the 
nuclear kinetic energy and is a Rieman- 

nian metric (Eubini-Study metric) defined as^^ 

Qrv = - |4>i?)(4>i?||Vzy4>i?). (32) 

This is directly analogous to the result already known in 
the Born-Oppenheimer approximation.The gra¬ 
dient of fgeo gives a geometric contribution to the elec¬ 
tric field acting on the nuclei.The metric g^y and 
Berry curvature B^y can be unified into a quantum geo¬ 
metric tensor.Switching to arbitary collective nuclear 
coordinates the quantum geometric tensor is 

= {Or^rII - \^R){^R\\dy^R). (33) 

The real part of T^y is the metric g^^y while the imaginary 
part is 1/h times the Berry curvature B^y. 

In the adiabatic approximation, the potential fgeo 
diverges with the inverse square distance from a con¬ 
ical intersection between two energy surfaces.By 
repelling the nuclei from the vicinity of the conical 
intersection where the adiabatic approximation breaks 
down, the potential fgeo enhances the accuracy of that 
approximation.^^ In line with what we found for the vec¬ 
tor potential and potential energy surface in previous sec¬ 
tions, the singularity of the adiabatic is smoothed 
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FIG. 8: Geometric contribution ^geo to the potential energy. 


out in the exact quantity. Figure 8 shows that the diver¬ 
gent adiabatic potential is rounded off to a smooth finite 
peak when the metric is evaluated with the exact 
electronic function instead of the Born-Oppenheimer 
function 


VII. CONCLUSIONS 

The topological Longuet-Higgins phase^ and accompa¬ 
nying vector potential^ are observable in numerous exper¬ 
iments, e.g. spectroscopy of triatomic molecules^^’^^’^^ 
and dynamical Jahn-Teller defects in bulk crystals^^”^^. 
The identification of such topological phases has always 
been based on the Born-Oppenheimer approximation, 
leaving open the possibility that their topological charac¬ 
ter is an artifact of that approximation which would not 
survive in an exact calculation. By identifying a specific 
case where a nontrivial Longuet-Higgins phase of tt does 
indeed vanish in an exact calculation based on the exact 
electron-nuclear factorization, recent work has amplified 
this uncertainty.^^ In the model of Ref. 14, the molecular 
geometric phase only takes the adiabatic value tt if there 
is a cusp (nonanalyticity) in the potential energy surface. 
Since the exact potential energy surface is smooth and a 
cusp is only recovered in the limit M ^ oc, the molec¬ 
ular geometric phase jumps discontinuously from tt to 0 
when the mass is decreased from infinity to a large finite 
value, i.e. when the nuclear kinetic energy is turned on. 

However, it would be puzzling if that behavior were to 
occur in the classical models of pseudorotating molecules 
because in those models the Longuet-Higgins phase is a 
topological invariant identifiable from qualitative global 
properties of the electronic Born-Oppenheimer wavefunc- 
tion 4>^^(r) far away from the conical intersection^^ and 
therefore robust to perturbations. For large M, the per¬ 
turbation induced by turning on the kinetic energy is 
localized near the conical intersection and should not af¬ 
fect its global properties. Therefore, the exact ^R{r) is 
almost everywhere similar to 4>^^(r) and the geometric 


phase should not be expected to jump from tt to 0. 

To resolve this discrepancy, we have applied the exact 
factorization to a model pseudorotating molecule which 
is closer in spirit to the original example of Herzberg 
and Longuet-Higgins.^ We have been able to answer the 
two questions raised in the introduction. First, when the 
Born-Oppenheimer factorization is replaced by the exact 
factorization, the conical intersection of the adiabatic po¬ 
tential energy surfaces is smoothed out, and second, the 
Longuet-Higgins phase becomes a path-dependent 1/(1) 
geometric phase. This is our main result: quantities that 
were discrete topological invariants in the adiabatic ap¬ 
proximation change into geometric quantities in the exact 
factorization. Since geometric phases can take any value 
between 0 and 27r, not only 0 or tt, the molecular geo¬ 
metric phase of the current-carrying state |T+) decreases 
continuously from tt to a value slightly less than tt as M 
is reduced from infinity to some large value. The molec¬ 
ular geometric phase of a current-carrying state remains 
finite even though the exact potential energy surface is 
everywhere smooth, proving that nonanalyticity is not 
a necessary condition for nonzero molecular geometric 
phase as previously believed. Whether there might be 
other topological contributions to the molecular geomet¬ 
ric phase, e.g. from nodes of x(R), is an open question. 

Unlike the Longuet-Higgins phase, which is only path 
dependent insofar as it depends on the winding number of 
the path around the conical intersection, the exact molec¬ 
ular geometric phase (like the Berry phase) is truly path 
dependent. The integral ^ • dH^ is the expression 

for the flux of a smeared-out fiux tube through a sur¬ 
face bounded by the path C. Since the surface does not 
enclose all of the flux of the smeared-out flux tube, the 
integral is not quantized to 0 or tt. In the limit M ^ oc, 
the smeared-out flux tube shrinks to a line and the exact 
geometric phase approaches the adiabatic value 0 or tt. 
Degeneracy is the crucial factor allowing us to construct 
current-carrying states with nontrivial Berry curvature. 
This is the essential difference with respect to the model 
studied in Ref. 14, which does not have a degeneracy fol¬ 
lowing from pseudorotational symmetry because two of 
the three nuclei are held fixed. 

The model studied here provides instructive examples 
of nontrivial molecular geometric phase and induced vec¬ 
tor potentials and exact 2D potential energy surfaces. 
Choosing a real-valued state from the degenerate ground 
state manifold leads to an exact potential energy sur¬ 
face that necessarily breaks the threefold symmetry of 
the model Hamiltonian. This is in contrast to the adi¬ 
abatic case, where there is a single symmetric potential 
energy surface for all states belonging to a degenerate 
manifold. On the other hand, threefold symmetry can be 
preserved by choosing a complex current-carrying ground 
state, but then the nuclear Schrodinger equation must 
contain a nontrivial vector potential. Thus, it appears 
that choosing a current-carrying state leads to a poten¬ 
tial energy surface that is closer to the familiar smooth 
and symmetric adiabatic potential energy surfaces. 
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Nontrivial induced vector potentials can occur in any 
molecular system due to the degeneracy associated with 
rotational symmetry, although in some cases they may be 
negligible. It will be necessary to account for induced vec¬ 
tor potentials in some time dependent problems, e.g. ex¬ 
citations to degenerate excited states with significant nu¬ 
clear currents. Methods to exploit the exact factorization 
in coupled electron-nuclear dynamics are under active 
development,^^ but so far have only been applied to one¬ 
dimensional systems, where is trivial. The physical 
role of the induced vector potential is to provide a contri¬ 
bution to the nuclear current which the gradient 

term filmy*is not able to provide. 

Note added in manuseript - Englman has applied 
the exact factorization to the model of Longuet-Higgins 
and coworkers^^ in a paper^^ that appeared after our 
manuscript was submitted. The model is a special case 
of the one studied here, but Ref. 76 considered only the 
adiabatic (Longuet-Higgins) phase, 0 or tt, and not the 
exact geometric phase defined by Eqs. (21) and (22). 
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The matrix elements of the internuclear repulsion obey 
the selection rules Am = 0 and AN = 0 , ± 1 . If AN = 0 , 
we have 


{nm\Vnn\nm) = —(n + 1 ). 


and if AN = ±1, we have 


{nim\Vnn\n2m) = N2){\m\ + max{Ni, N2)). 


The electron-nuclear coupling H^n is the only term 
which couples states with different values of m. To sim¬ 
plify the evaluation of its matrix elements, we first change 
to the nuclear coordinate representation 

n27r ^ nOO 

{ainimi\Hen\a 2 n 2 m 2 ) = — QdQp niTTii (Q) 

Jo Jo 

(Al) 


Appendix A: Hamiltonian matrix elements 


All radial integrals can be evaluated analytically with 
the recursion relations for the Laguerre polynomials.^^ 
We start with the first term of Eq. ( 11 ), which obeys 

the selection rule Am = m2 — mi = 0 . We distinguish 
three cases: (i) AN = N2 — Ni = 0 , (ii) AN = ±1 and 
(hi) |AA| > 1, where Ni = (n^ — |m|)/2. In case (i) we 
find 

(nm|T^^^|nm) = hft(^N ^ 

in case (ii). 


(nim|T^^^|n2m) 
and in case (hi). 


hfl (max(A/’i, ^2))^/^ 

2 ^y\m\ +max(Ai, A2)’ 




{nim\T'^^'^\n2m) = — —\m\ 


/max(Ai,A2)! 
min(A/’i, A2)! 



min(Ai, A2))! 


(|m| +max(Ai, A2))!* 


We have defined the frequency Q = The matrix 

elements of the second term of Eq. ( 11 ), obey the 

selection rule Am = 0 . They are the same as those for 
case (hi) of but with opposite sign, 

(nim|T^^^|n2m) = —(nim|T^^^|n2m) [case (hi)]. 


Using Eqs. (6) and ( 7 ), = {aiQr]\Hen\ci2Qr]) 

is divided into the following three contributions: 


H^^{Qr]) = I ^ J gQcosf] + K ^Qsinr^. (A 2 ) 


/, J and K are 8x8 matrices reported in Appendix B. 
The matrix elements of Qcosr] and Q sin 77 can be re¬ 
covered from the following matrix elements of Qe'^^ and 
which were derived in Ref. 17 : 


(n, mjQe -h 1, m + 1) = (n + 1, m + IjQe^^ln, m) 


n + m + 2 


(n, m\Qe — 1, m + 1) = (n — 1, m + IjQe^^ln, m) 



Appendix B: Electron-nuclear coupling matrices 


As a function of (Q,?^), the electron-nuclear coupling 
can be represented as shown in Eq. (A2). To see how 
parity symmetry affects the matrices /, J and A, it is 
convenient to express them in the following real-valued 
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electronic basis: 


and 


W) 

W) 

\c') 

\d') 

|e') 

I/') 

|5') 

\h’) 


^i\a) + \b)) 
^(|c) + M)) 

;i(|e)-|/)) 

;^(|e) + l/)) 

-^15) 

-\h). 


(Bl) 


The states |a'), |c'), |e') and \g') are odd with respect to 
the reflection Qs ^ — Qs, while |6'), jd'), |/') and |h') 
are even. In the ja') basis, we have 

I = diag(—3, —3, +3, +3,0,0,0,0), 


/ 0 

-1 

0 

0 

0 

1 

2 

0 

+1 

\ 

-1 

0 

0 

0 

1 

2 

0 

x /3 

2 

0 


0 

0 

0 

+1 

0 


0 



0 

0 

+1 

0 


0 

1 \/3 
^ 2 

0 


0 

1 

2 

0 


0 

0 

0 

2 


1 

2 

0 


0 

0 

0 

0 

0 


0 

2 

0 

1 \/3 

2 

0 

0 

0 

0 



0 

+1 

0 

2 

0 

0 

0 

/ 


/ 1 0 0 0 +1 0 0 \ 

0 -1 0 0 0 -i 0 -i 

0 0 -1 0 0 0 

0 0 0 1 0 0 -i 

+ i 0 0 0 0 0 0 

0 0 +1 0 0 0 -2 
0 00000 
V 0 -i 0 -i 0 -2 0 0 / 


The matrix J couples states with the same parity because 
cos 7^ is even under reflection. The matrix K couples 
states of different parity because sinr^ is odd. 
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